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Abstract. The accretion column in magnetic Cataclysmic Variables may have a not negligible Thomson optical 
depth. A fraction of the thermal radiation from the post-shock region may therefore be scattered - and then 
polarized - before escaping the column. Moreover, part of the thermal radiation is reflected - and again polarized 
- by the White Dwarf surface. We show that X-ray polarimetry can provide valuable, and probably unique 
informations on the geometry and physical parameters of the accretion column by calculating, by means of 
Monte Carlo simulations, the expected polarization properties of magnetic CVs as a function of the geometrical 
parameters (assuming a cylindrical geometry) and the Thomson optical depth of the column. We find that degrees 
of polarization as high as about 4% can be present, and apply our calculations to the archetypal magnetic CV, 
AM Herculis. 
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1. Introduction 

The magnetic field in some subclasses of Cataclysmic 
Variables (notably in Polars and in at least a fraction 
of intermediate Polars; see Warner 1995 for a complete 
overview on CVs) is strong enough to channel the accret- 
ing matter along the field lines. For a dipolar field, this 
means that the accretion occurs via accreting columns on 
one or both the magnetic poles; the disalignement between 
the spin and magnetic axes results in pulsed emission. 

Hard X-rays are then produced by optically thin ther- 
mal line and continuum emission in the so-called post- 
shock region, where temperatures can reach values as large 
as several tens of keV (e.g. Frank et al. 1992; Cropper et 
al. 1999). Thermal emission is expected to be unpolarized 
(see next section); however, the Thomson optical depth in 
the accretion column, while probably less than unity, may 
be not negligible. Indeed, Hellier et al. (1998) found signif- 
icant broadening in the iron Ka line of several magnetic 
CVs, which they interpreted as due to Compton broad- 
ening. Thomson scattered radiation is polarized, provided 
that the geometry is not spherical. It is therefore to be 
expected that the hard X-ray emission in magnetic CVs 
is polarized, with the net polarization degree increasing 
with the Thomson optical depth of the accretion column. 

In this paper we calculate, by means of Monte Carlo 
simulations, the polarization properties of the accretion 
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column in magnetic CVs. In Sec. 2 the numerical code 
is described, while the results (including an application 
to the archetypal Polar, AM Herculis) are presented in 
Sec. 3. Results are then summarized, and observational 
perspectives discussed, in Sec. 4. 



2. Calculations 

The Monte Carlo technique adopted in our code is de- 
scribed in Matt (1993) and Matt et al. (1989; 1996). A 
full Compton scattering matrix is adopted (McMaster 
1961). However, in the energy range under consideration 
( ^ 20 keV) , the scattering matrix basically reduces to the 
Rayleigh one. We assumed that primary photons (i.e. pho- 
tons before scattering) are unpolarized. Bremsstrahlung 
photons have a polarization, i.e. electric, vector perpen- 
dicular to the interaction plane, and the random velocity 
field ensures that the radiation has a null net polarization. 
In practice, in the code each photon is assigned a random 
polarization vector. After emitted, the photon path is fol- 
lowed until it leaves the accretion column. We assumed 
total ionization of the matter, i.e. no photoelectric ab- 
sorption. 

We assumed, for simplicity, a cylindrical geometry, 
constant density (and therefore emissivity) and zero tem- 
perature of the electrons. While the first assumption, as 
any geometrical assumption, is critical for polarization 
calculations, the second assumption is not very relevant. 
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Fig. 1. The degree of polarization, P, as a function of fi = cos6 for four values of H/rc- In each panel, five values of 
tt are shown: 0.05, 0.1, 0.3, 0.5 and 1 (in increasing order of the absolute value of P). 



the important parameter being the total scattering optical 
depth. The third assumption is apparently very strong - 
typical maximum temperatures of magnetic CVs arc 10- 
20 keV for polars, and up to 30-40 keV for Intermediate 
Polars. However, especially for polars and for unpolarized 
primary emission, the zero temperature approximation is 
still reasonably good as far as polarization properties (see 
e.g. Poutanen & Vilhu 1993) and intensity distributions 
are concerned. 



We calculated the polarization properties as a func- 
tion of two parameters: the radial Thomson optical depth 
(tt), and the H/rc ratio, where rc is the radius and H 
the height of the cylinder. Simple modeling of the accre- 



tion column gives, for these two parameters, the following 
equations (e.g. Prank et al. 1993): 

H/rc ^ ^-^M^^'fhML^iR'X (1) 

TT Ommefl^M-l.RZl, (2) 

where Miq is the accretion rate in units of 10^^ g s~^, 
/_2 is the fractional area of the column on the White 
Dwarf surface in units of 10~^ (/=^c/4-R4d)' ^wd,i the 
mass of the White Dwarf in units of the solar mass, and 
Rwd,9 the radius of the White Dwarf in units of 10^ cm. 
It is worth noticing that tt increases with the accretion 
rate, and therefore we expect a larger polarization degree 
when the sources are in their high states. It must also be 
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Fig. 2. The same as in the previous figure, for other four values of H/rc- 



recalled that eq.^is valid for bremsstrahlung cooling only. 
If cyclotron cooling is not negligible, eq. ^ gives only an 
upper limit to H/rc (Frank et al. 1993).. 

To decide the range of parameters to be explored, we 
assumed i?„ci, 9=1 and Mwd.i=^i and took the results on 
the accretion rate and column radius given in Cropper 
et al. (1998). These parameters were obtained by fitting 
GINGA spectra with a more complex model, and therefore 
our procedure is not fully self-consistent, but certainly 
good enough for an order-of-magnitude estimate. 

tt results to be typically in the range 0.05-1. This 
is therefore the range we will explore (in any case, for 
Tt ^ 0.05 the polarization degree is neghgible). 

H/rc is typically in the range 1-10 (apart for sources 
in low states, when it can became much larger but the 
optical depth is very low). Given the fact that H/rc is an 



upper limit, we extend the explored range down to 0.1. It 
must be noted that for H/rc larger than 10 and smaller 
than 0.1, the polarization properties will not change much 
any further because the geometry 'saturates', apart from 
very small (in the former case) or large (in the latter case) 
inclination angles when border effects are still present. 

Finally, it is important to note that the two parameters 
are not independent of each other, but H/rc increases 
with / and decreases with M, while tt has the opposite 
behaviour. We therefore expect large radial optical depths 
for flattened geometrical configurations. For the sake of 
completeness, however, we will explore the full parameter 
space. 

1 billion input photons are used for each simulation. 
The emerging photons are stored in angular bins, with 
A/i=0.1 (/i — cosO, where 9 is the inclination angle mea- 
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Fig. 3. The flux (in arbitrary units) as a function of /i = cos9 for four values of H/rc- In each panel, five values of tt 
are shown: 0.05, 0.1, 0.3, 0.5 and 1 (in increasing order of anisotropy). 



sured from the cylinder axis, i.e. the angle bewteen the 
magnetic axis and the line-of-sight). 



3. Results 

The hard (i.e. above ~1 keV) X-ray emission of mag- 
netic CVs consists not only of plasma (bremsstrahlung) 
emission but also of a reflection component, arising from 
the illumination of the WD surface (e.g. van Teeseling et 
al 1994; Beardmore et al 1995; Matt 1999 and references 
therein) . 

Let us first present the polarization properties of the 
accretion column alone, and afterwards we will discuss the 
contribution of the reflection component. 



3.1. Emission from tlie accretion column 

As long as the scattering matrix docs not differ too much 
from the Rayleigh one, the polarization properties of the 
X-ray radiation from the accretion column are energy- 
independent. For obvious symmetry reasons, the polariza- 
tion is expected to be either parallel or perpendicular to 
the projection of the cylinder axis onto the plane of the 
sky. In the former case, as customary when dealing with 
axisymmetric geometries, we will conventionally assume 
that the polarization degree, P, is negative, while positive 
of course in the latter case. 

In Figs. Handle the polarization degree as a function 
of /i for different values of tt and H/rc, is shown. For 
H/rc ^ 1, the polarization is negative, as expected for a 
disc-like geometry and tt less than a few (e.g. Sunyaev 
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Fig. 4. The same as in the previous figure, for four other values of H/rc- 



& Titarchuk 1985). Of course, P increases (in absolute 
value) with tt, as more and more photons are scattered 
(and therefore polarized) before escaping. P at first de- 
creases with increasing H /rc, becoming very small when 
this parameter is equal to 1, not suprisingly as in this 
case the geometry is rather symmetric. Further increas- 
ing Hjrci the polarization degree (now positive) starts to 
increase again. 

In Figs. inland^ the flux (in arbitrary units) as a func- 
tion of /i is instead shown. The flux anisotropy of course 
increases with increasing tq- 

It is worth recalling that the X-ray polarization is ex- 
pected to be phase-dependent, because of the dependence 
on the spin phase of 0, which is basically the angle between 
the line of sight and the magnetic dipole axis. It must be 
also remarked that 9 can become larger than 90°, in which 



case the accretion column is partly occulted by the WD 
surface. The degree of polarization in this case depends 
of course on the level of occultation, but it is expected 
to be larger than in absence of occultation, because the 
geometry becomes less symmetric. 



3.2. The reflection component 



The reflection component is polarized parallel to the pro- 
jected cyhnder axis (Matt et. 1989; Matt 1993), i.e. it is 
negative. The degree of polarization as a function of 9 is 
shown in Fig. |31 again obtained by means of Monte Carlo 
simulations (Matt et al. 1989; photons are integrated over 
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Fig. 5. The polarization degree of the reflected compo- 
nent, as a function of /i. 

Fig. 6. The intensity (upper panel) and degree of polar- 
thc 1-20 keV energy range). The net polarization is then ization (lower panel) expected for AM Herculis, calculated 



given by: 



Ppr{lJL)+r{E,ll)Pref{E,tl) 



(3) 



where Ppr and Pref sue the polarization degrees of the 
primary and reflected emission, respectively (taken with 
their signs), and r(E, ii) is the ratio between the primary 
and reflected fluxes (see e.g. Matt 1999). 

3.3. The expected X-ray polarization of AM Herculis 

As an application of our calculations, let us discuss 
the expected phase-dependent X-ray polarization of the 
archetypal polar, AM Herculis. Assuming the system pa- 
rameters given by Cropper (1988), i.e. P—61° and i=30° 
(where /3 is the angle between the magnetic and spin axes, 
and i the angle between the spin axis and the line-of- 
sight), 9 never exceeds significantly 90°, so our calcula- 
tions can be safely applied to the whole spin period. 

The intensity and degree of polarization of the accre- 
tion column on the main magnetic pole, without the reflec- 
tion component (and therefore valid for energies less than 
a few keV), as a function of the spin phase, </>, are shown in 
Figs. Eland In Fig.El H/rc=lO and rT=0.05 have been 
assumed, according to the value derived from the Cropper 
et al. (1998) estimates from ASCA data. However, during 
the ASCA observation the source was in an intermediate 
X-ray state; larger optical depths, and lower H/rc, are 
expected at higher luminosities. Moreover, it must be re- 
called again that the H/rc values estimated from eq.^are 



without the reflection component, as a function of the spin 
phase. H/rc—W and tt=Q.05. 
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Fig. 7. The intensity (upper panel) and degree of polar- 
ization (lower panel) expected for AM Herculis, calculated 
without the reflection component, as a function of the spin 
phase. H/rc— 0.5 and tt=0.5. 
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upper limit, as cyclotron cooling is not taken into account. 
We therefore calculated the phase-dependent polarization 
properties also for the (more favourable) case of H/rc—0.5 
and tt=0.5 (Fig.[7|. 

In the first case (Fig. EJ, the polarization degree is 
low (less than 1%) because of the low optical depth. The 
flux is basically phase-independent, differently from what 
found in intermediate states by BeppoSAX and ASCA 
(Ishida et al. 1997, Matt et al. 2000). Unless the real op- 
tical depth is larger than derived from eq. 2, this means 
that either the assumed geometrical and physical parame- 
ters of the accreting column are oversimplified (it must be 
recalled that we have assumed a single temperature and 
density along the column), or that there is substantial 
phase-dependent absorption. Studying deviations of the 
measured polarization from the expected one may help 
distinguishing between the different hypotheses. 

The phase-dependent flux shown in Fig. [T] is instead 
roughly similar to that observed by BeppoSAX during 
a high state of the source (Matt et al. 2000); the ob- 
served amplitude of the modulation is however larger then 
predicted, again suggesting a more complex gemetrical 
and/or physical situation than assumed in our calcula- 
tions. The polarization degree is expected to be as large 
as 1.5%, increasing when the flux decreases. 

The phase-dependent degree of polarization, includ- 
ing the reflection component as estimated by Matt et al. 
(2000) during a high state of the source, is shown in Fig.|Sl 
for two energy bins (5-10 keV and 10-15 keV); Rc/H=10 
and rT=0.05 (upper panel) and Rc / H=0.5 and tt=0.5 
(lower panel). 
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Fig. 8. The degree of polarization, including the reflec- 
tion component, expected for AM Hcrculis as a function 
of the spin phase, for two energy bins (5-10 keV, filled 
circles, and 10-15 keV, crosses). Upper panel: H/rc=lO 
and rT=0.05; lower panel: H/rc=0-5 and tt=0.5. 



4. Summary 

We have calculated the polarization properties of the ac- 
cretion column in magnetic CVs. Polarization arises be- 
cause a fraction of the thermal radiation can be Thomson 
scattered before escaping the column. The polarization 
degree can be as high as ^4%, and depends on the angle 
between the magnetic field and the line-of-sight (which 
of course varies with the spin phase), which is often 
well known (e.g. Cropper 1988). It depends also on the 
Thomson optical depth, tt, and on the ratio between the 
radius and height of the accretion column. If one of these 
two parameters can be independently estimated (e.g. tt 
from the iron line broadening, e.g. Hellier et al. 1998), 
than the other can in principle be deduced by polariza- 
tion measurements. In particular, as the polarization de- 
gree is negative (positive) for H /rc less (greater) than one, 
while the polarization of the reflection component (whose 
importance increases with energy) is always negative, an 
increase (decrease) of the polarization degree with energy 
is expected in the former (latter) case. 

In this paper we have discussed the continuum emis- 
sion. However, iron lines provide a significant fraction of 
the total X-ray flux. Recombination lines emitted in the 
accretion column may suffer not only Compton scatter- 
ing but also resonant scattering, so the effective optical 



depth in the line is larger than in the continuum. Line 
photons are therefore likely to be more polarized than con- 
tinuum photons. On the contrary, the fluorescent neutral 
line emitted by the White Dwarf surface is unpolarized, 
at least in the line core. 

Degrees of polarization of the order of one percent are 
within the detection capabilities of the new generation X- 
ray polarimeters based on the photoelectric effect (Costa 
et al. 2001) when coupled with large enough X-ray tele- 
scopes (Costa et al. 2003), at least for the brightest mag- 
netic CVs in high state (when their flux can be as high 
as several millicrabs; for instance the 2-10 keV flux of 
AM Herculis in high state is about 10~^° erg cm~^ s~^, 
Matt et al. 2000). Moreover, the spectral resolution should 
be good enough to search for different polarization degrees 
in the emission iron lines. Bright magnetic CVs should 
therefore be added to the traditional lists of targets for 
future polarimetric missions. 
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